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Abstract. We consider the weakly asymmetric exclusion process on a bound- 
ed interval with particles reservoirs at the endpoints. The hydrodynamic limit 
for the empirical density, obtained in the diffusive scaling, is given by the vis- 
cous Burgers equation with Dirichlet boundary conditions. In the case in which 
the bulk asymmetry is in the same direction as the drift due to the boundary 
reservoirs, we prove that the quasi-potential can be expressed in terms of the 
solution to a one-dimensional boundary value problem which has been intro- 
duced by Enaud and Derrida 1151 . We consider the strong asymmetric limit of 
the quasi-potential and recover the functional derived by Derrida, Lebowitz, 
and Speer |14l for the asymmetric exclusion process. 
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1. Introduction 

The study of steady states of non-equilibrium systems has motivated a lot of 
works over the last decades. It is now well established that the steady states of 
non-equilibrium systems exhibit in general long-range correlations and that the 
thermodynamic functionals, such as the free energy, are not local nor additive. 

The analysis of the large deviations asymptotics of stochastic lattice gases with 
particle reservoirs at the boundary has proven itself to be an important step in 
the physical description of nonequilibrium stationary states and a rich source of 
mathematical problems. We refer to [51 [T3| for two recent reviews on this topic. 

We consider a boundary driven one-dimensional lattice gas whose dynamics can 
be informally described as follows. Fix an integer iV > 1, an external force E in 
M and boundary densities < p_ < p+ < 1. At any given time each site of the 
interval {—N -I- 1, . . . , iV — 1} is cither empty or occupied by one particle. In the 
bulk, each particle attempts to jump to the right at rate 1 + E/2N and to the left at 
rate 1 — E/2N. To respect the exclusion rule, the particle jumps only if the target 
site is empty, otherwise nothing happens. At the boundary sites ±(iV— 1) particles 
are created and removed for the local density to be p±: at rate p± a particle is 
created at ±(A^ — 1) if the site is empty and at rate I — p± the particle at ±{N — 1) 
is removed if the site is occupied. 

The dynamics just described defines an irreducible Markov process on a finite 
state space which has a unique stationary state denoted by fi^. Let (p± := 
log[p±/(l — P±)] be the chemical potential of the boundary reservoirs and set 
Eq := — (y9_)/2. When E = Eq, the drift caused by the external field E matches 
the drift due to the boundary reservoirs, and the process becomes reversible. 

In the limit | oo, the typical density profile p^ under the stationary state p^ 
can be described as follows. For each E < Eq there exists a unique Je < such 
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that 

1 /■"+ . 1 



dr ——- — = 1 



2 A_ Exir) - .1 



E 



where x is the mobiUty of the system: xi'^) ~ cl{1 — a). The profile is then 
obtained by solving 

p'e -Ex{Pe) = - Je 
with the boundary condition 1) = p_. 

In the same limit N ] oo, the probability of observing a density profile 7 different 
from Pe can be expressed as 

{7} exp{-Wi,(7)} • (1-1) 

The large deviations functional Ve, which also depends on /9_, p+, is an extension 
of the notion of free energy to the context of non-equilibrium systems. 

The free energy of a boundary driven lattice gas has first been derived for the 
symmetric simple exclusion process by Derrida, Lebowitz and Speer [14j based 
on the so called matrix method, introduced by Derrida, which permits to express 
the stationary state p,^ as a product of matrices. Bertini et al. [3] derived the 
same result through a dynamical approach which we extend here to the weakly 
asymmetric case. 

We consider only the situation E < Eq for the bulk asymmetry to be in the 
same direction as the drift due to the boundary. The reversible case E = Eq lacks 
interest because the stationary state is product and does not exhibit long range 
correlations. In contrast, the analysis of the quasi-potential Ve for E > Eq, not 
treated here, appears a most interesting problem. For instance, a representation of 
Ve as a supremum of trial functionals analogous to p.l4p below seems to be ruled 
out. 

In the boundary driven weakly asymmetric exclusion process, for E < Eq, the 
quasi-potential takes the following form: 

Ve{j) ■■= j rfii|7log7+(l-7)log(l-7) + (l-7)</j-log(l + e^) 



-1 

1 

+ E 



if' \og<p' - - E)\og[>p' - E)] - Ae] , (1.2) 



where A^; is the constant given by 



7+ 



Ae := log(-j£;) + \ I dr ^ log 
2 Ex{r) 



Exir) 



Je 

and where ip is the unique solution of the Euler-Lagrange equation 

ip" 1 _ 

ip'iip' -E) ^ l + e'P ~ ^ 

satisfying v?(±l) = p±, ip' > max{0, i?}. 

This result, stated in a different form, has been proved by Enaud and Derrida 
[TS] based on the matrix method. We prove this result in Section 2] below by the 
dynamical approach introduced in [3] 2] . We also show that the quasi-potential is 
convex and lower semi- continuous. 

In section [51 we show that Ve F-converges, as ii^ J, —00, to the free energy 
of the boundary driven asymmetric exclusion process, first derived by Derrida, 
Lebowitz and Speer [13] . This asymptotic behavior is somewhat surprising since the 
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hydrodynamic time scales at which the weakly asymmetric exclusion process and 
the asymmetric exclusion process evolve are different. We also prove convergence 
of the solutions of the Euler-Lagrangc equations as the external force E diverges. 

The dynamical approach followed here permits to compute the fluctuation prob- 
abilities in great generality, in any dimension and for a large class of processes. 
However, it is only in dimension one and for very few interacting particle systems 
that an explicit expression of type (|1.2p is available for the non- equilibrium free 
energy Ve- 

2. Notation and Results 

The boundary driven weakly asymmetric exclusion process. Fix an integer 
TV > 1, S e M, < p_ < p+ < 1 and let Aat {-N + 1, . . . , TV - 1}. The 
configuration space is Sat := {0, 1}"^™; elements of S^r are denoted by 77 so that 
r]{x) = 1, rcsp. 0, if site x is occupied, resp. empty, for the configuration rj. We 
denote by a^'^rj the configuration obtained from rj by exchanging the occupation 
variables 7/(2;) and ri{y), i.e. 

(viy) z = x 
ri{x) \i z = y 
ri{z) if z ^ x,y, 

and by a^ij the configuration obtained from 77 by flipping the configuration at x, 
i.e. 

. J 1 - iri{x) if z = x 

I rj{z) II z ^ X. 

The one-dimensional boundary driven weakly asymmetric exclusion process is 
the Markov process on S ^ whose generator L n can be decomposed as 

Ln = Lq^n + L^^N + i+^jv , (2.1) 
where the generators Lq^n, L-^n, L+^n act on functions / : Sat R as 

{Lo,Nf){ri) ^ -Y E e-^/(2^)["(-+i)-''(-)l[/(a-'-+ir;)-/(r;)] , 

x = -N+l 

(L_,A./)(7y) = _c_(r;(-iV + l))[/(a-^+i^)-/(^)] 

{L+Mfm - — c+(7/(7V ~ 1)) [/(a^-^ry) - /(ry)] 
where c± : {0, 1} — > M are given by 

c±(C) p±eT^/(2^'(l-C) + (l-p±)e±^/('^)C. 

Notice that the (weak) external field is {2N) and, in view of the diffusive scaling 
limit, the generator has been speeded up by N"^ . We denote by 774 the Markov 
process on S^r with generator Ln and by its distribution if the initial configu- 
ration is rj. Note that is a probability measure on the path space D(M-|-, S^y), 
which we consider endowed with the Skorohod topology and the corresponding 
Borel (T-algebra. Expectation with respect to P^ is denoted by E^. 

Since the Markov process -qt is irreducible, for each > 1, iiJ G R, and < p_ < 
p+ <1 there exists a unique invariant measure /x^ in which we drop the dependence 
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on p± from the notation. Let if± := log[p±/{l — p±)] be the chemical potential of 
the boundary reservoirs and set Eq '■= — ip-)/2. A simple computation shows 
that if _E = i^o then the process rjt is reversible with respect to the product measure 

n 7^;^ (2-2) S 

where 

_N , ^ iV - X N + X 

On the other hand, for E ^ Eq the invariant measure /i;^ cannot be written in a 
simple form. 

The dynamical large deviation principle. We denote by u G [1,1] the macro- 
scopic space coordinate and by (•, •) the inner product in L2([— 1, 1], du). We set 

M:= {peLoo{[-l,l\,du) : 0<p<\} (2.3) | dcm 

which we equip with the topology induced by the weak convergence of measures, 
namely a sequence {p"} C M converges to p in if and only if {p",G) — > {p, G) 
for any continuous function G : [— 1, 1] ^ M. Note that is a compact Polish space 
that we consider endowed with the corresponding Borel cr-algebra. The empirical 
density of the configuration 77 G Sat is defined as 7r^(r/) where the map tt^ : E^v 
M is given by 

N-l . 

n^{v){u):= E ,ix)l{[^-—,^ + —)}{u), (2.4) ^ 

in which 1{j4} stands for the indicator function of the set A. Let {r/^} be a sequence 
of configurations with ry^ G Sat. If the sequence {TT^{ri^)} C Ai converges to p 
in as — > cxD, we say that is associated the macroscopic density profile 

p e M. 

Given T > 0, we denote by D{[0,T];M) the Skorohod space of paths from 
[0,T] to M equipped with its Borel a-algebra. Elements of D[[0,T], M) wih be 
denoted by tt = Trt{u) and sometimes by ■n{t,u). Note that the evaluation map 
-D([0, r];A^) 9 TT i-^ TTt G is not continuous for t G (0, T) but is continuous 
for t = 0,T. We denote by tt^ also the map from /^([O, T]; Sjv) to D([0,r];A4) 
defined by -k^ {r}.)t := TT^{T]t). The notation TT^{t,u) is also used. 

Fix a profile 7 G and consider a sequence {77^ : -/V > 1} associated to 7. 
Let T]^ be the boundary driven weakly asymmetric exclusion process starting from 
7]^. In [12l[m[20] it is proven that as ^ 00 the sequence of random variables 
{n^iv^)}, which take values in V{[0,T];M), converges in probability to the path 
p = pt{u), {t,u) G [0,r] X [—1, 1] which solves the viscous Burgers equation with 
Dirichlet boundary conditions at ±1, i.e. 

dtp + I Vx(p) - ^ Ap 

Pti±i) = p± (2-5) hq^i 

.Po{u) =^ 7(w) 

where x '■ [0, 1] K+ is the mobility of the system, x(a) = a(l — a), and V, resp. 
A, denotes the derivative, resp. the second derivative, with respect to u. In fact the 
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proof presented in [121 [18] is in real line, while the one in [20] is on the torus. The 
arguments however can be adapted to the boundary driven case, see [Ml [171 [19] for 
the hydrodynamic limit of different boundary driven models. 

A large deviation principle for the empirical density can also be proven following 
[20l [22l [23] , adapted to the open boundary context in ^ . In order to state this 
result some more notation is required. Fix T > and let fiy = (0,T) x (—1,1), 
r^T = [0, T] X [—1, 1]. For positive integers m, n, we denote by C™"'"{^It) the space 
of functions G = Gt{u) : Qt ^ K with m derivatives in time, n derivatives in space 
which are continuous up the the boundary. We improperly denote by C^'"(f2T) 
the subset of C™'"(r2T) of the functions which vanish at the endpoints of [—1, 1], 
i.e. G e C™^"(fV) belongs to C™'"(f^7) if and only if Gt(±l) = 0,te [0,T]. 

Let the energy Q : D{[0,T],A4) [0, oo] be given by 

QW = 

sup < dt 



H 



duTT{t,u){VH){t,u) - - / dt [ duH{t,uf x{T^{t,u))\ 
2 Jq J_i J 



where the supremum is carried over all smooth functions H : Qt — * K with compact 
support. If Q(7r) is finite, tt has a generalized space derivative, Vtt, and 



Q{tt) ^ IJ^ dt J du 



Fix a function 7 G which corresponds to the initial profile. For each H in 
Cg^{Th), let Jh = JT,H.-y - D{[0,T],M) — > M be the functional given by 

Jh{7t) := {ttt,Ht)^{7,Ho)- f dt{TTt,dtHt) 

Jq 

~ \ I dt{7rt,AHt) + ^ [ dtS/Htil) - ^ f dtVHt{-l) 
^ Jo ^ Jo ^ Jo 

~ ^ [di{x{^tlVH,) - \j\t{x{^,),{VHtf). 
Let /t( • I7) : ^([O, T],M) — > [0, +00] be the functional defined by 



/y(7r|7) sup Jh{t^) ■ (2.6) | ±01 

The rate functional It{-\i) ■ D{[0,T],M) [0, 00] is given by 

/^(^|7) = I ^^W^) if SW <oo, (2.7) 

1^ 00 otherwise. ' ' 

Is is proved in [7], for any E in R, that the functional /t(-|7) is lower semicontinous, 
has compact level sets and that a large deviations principle for the empirical measure 
holds. 

Theorem 2.1. Fix T > and an initial profile ^ in Ai. Consider a sequence {rj^ : 
^ 1} of configurations associated to 7. Then, the sequence of probability measures 
{P^N ° {^^)~^ : N > 1} on D{[0,T]; M) satisfies a large deviation principle with 
N and good rate function /t(-|7). Namely, It{'\j) '■ -D([0, T]; A^) [0, 00] 
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has compact level sets and for each closed set C C D{[0,T]; M) and each open set 
O C D{\0,T]-M) 

IS llogP^„(^^ eC) < -inf /T(^|7) 

N—*oo iV ttGC 

lim llogP^„(^^ eO) > - inf /T(7r|7) . 

The quasi-potential. From now on we consider only the case E < Eq = (95+ — 
ip-)/2, where ip± = log[p±/(l — p±)]. Simple computations, which are omitted, 
show that the unique stationary solution G of the hydrodynamic equation 
(|2.5p can be described as follows. For each E < Eq there exists a unique Jb < 
such that 

1 rP+ 1 

- / dr — = 1 . (2.8) rciil 

2jp_ Ex{r)-JE ^ ^ ' ' 

The profile is then obtained by solving 

p'e-Ex{Pe)^-Je (2.9) H 

with the boundary condition P£;(— 1) = P-- Note that J_b/2 is the current main- 
tained by the stationary profile p^. The solution to (|2.9|) can easily be written in 
an explicit form, see [E]. We shall however only use, as can be easily checked, that 
is strictly increasing and that the inequality Je/ E > niaxrg[p_^p^] xi''') holds for 
E <0. 

Given E < Eq, the quasi-potential for the rate function It is the functional 
Ve ■ M ^ [0, +00] defined by 

VEip) := mf^ inf {/T(7r|ps) , ir £ D{[0,T];M) : ttt = p} (2.10) ^ 

so that Ve(p) measures the minimal cost to produce the profile p starting from p^. 

Recall that p^ is the unique invariant measure of the boundary driven weakly 
asymmetric exclusion process. The following result, which states that the quasi- 
potential gives the rate function of the empirical density when particles are dis- 
tributed according to p^ is proven in [S] in the case E = 0. Thanks to Theorem l2.11 
the proof applies also to the weakly asymmetric case. 

Theorem 2.2. For each E in M, the sequence of probability measures on A4 given 
by {Pe °{^^)^^} satisfies a large deviation principle with speed N and rate function 
Ve- Namely, for each closed set C C A4 and each open set O C M, 

llog^f(^^GC) < -infVEip) 

fim llogp^in^ eO) > -M VEip). 

In this paper we prove that the quasi-potential Ve can be expressed in terms of 
the solution to a one-dimensional boundary value problem. This result has been 
obtained in |15| by analyzing directly the invariant measure p^ through combi- 
natorial techniques; while we here follow instead the dynamic approach [3l [4] by 
characterizing the optimal path for the variational problem (j2.10p . 

For E < Eq, let C^+^([— 1, 1]) be the set of continuously differentiablc functions 
on [—1,1] with Lipshitz derivative and set 

^f;:-{^GCi+i([-l,l]) : (p(±l)=^±, (^'>0V£;} , (2.11) 
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Note that Se is a positive functional because a simple computation relying on (|2.9p 
shows that 



Se{p) > Qe{p,Ve) 



Grf 



where, given a,b ^ K, the notation aW h, resp. a Ab, stands for max{a,6}, resp. 
min{a,5}. Note that Te = Te' for E,E' < 0. 

For E < Eo, E ^ 0, let Ge ■■ M X Te ^ be given by 

Ge{p,^) ■■= J du{plogp+{l- p)log{l- p) + {l- p)^-log{l + e'P) 

+ ^[^' log ^' ^{^' -E) log{^' -E)\-Ae} , (2.12) 

where, by convention, OlogO = and Ae is the constant given by 

I rP+ 1 r EY(r)-i 

Ae log(-J£;) + - / dr-— -log 1-^1 . (2.13) 

2Jp_ Ex{r) I Je J I 1 

The right hand side is well defined because Je < and Je/E > maxrg[p_.p^] x('') 
for E <0. 

For E ~ Q, Qo : A4 X Te ^ M is defined by continuity as 
Qq{p,^) = 

J dM{plogp+(l-p)log(l-p) + (l-p)(/.-log(l + e^)+ log^' + l-^o} , 

where = log[(p+ - p_)/2] + 1. 

For E < Eq, define the functional S'^ : ^ R by 

Se{p) ■■= sup Qe{p,^p) ■ (2.14) I ssg 



/I ^ 
du|plog^ + (l-p)log-— ^1 (2.15) Vai^ 
\ ^ Pe ^- Pe^ ' ' 



f 07 



if ipE — log[P£;/(l - Pe)\- 

In the special case E = Eq, as already observed, the weakly asymmetric exclu- 
sion process is reversible and the stationary state Pe is a product measure. In 
particular, the rate functional Seq of the static large deviations principle for the 
empirical density can be explicitly computed. It is given by 

/I , 
d„|plog^ + (l-p)log-^| . (2.16) 
-1 Peo 1 - Pe„ J 

The Euler-Lagrange equation associated to the variational problem (j2.14p is 

7(^ + TT^ = ^- ^'-''^ S 

A function ip g J-e solves the above equation when it is satisfied Lebesgue a.e. 
Recalling that the stationary profile satisfies p.9p and p.Sp , it is easy to check 
that if p = then (pE solves p.l7p and Ge{pei ^e) ~ 0. 

The analysis of the quasi-potential for the boundary driven symmetric exclusion 
process, i.e. the case = of the current setting, has been considered in [i]. In 
particular it is there shown that Vq coincides with 5*0. We prove in this article an 
analogous statement for any E < Eq. 

t:S=V Theorem 2.3. Let E < Eq and Ve, Se : ^ [0, +oo] be the Junctionals defined 
m (|2T0ll . ([2T41) and ^JQ. 
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(i) The functional Se is bounded, convex, and lower semicontinuous on M. 

(ii) Fix E < Eq. For each p ^ A4 there exists in Te a unique solution to p.l7|) 
denoted by ^{p). Moreover 

Se{p) = max ^^(p, f) = GEip, Hp)) ■ (2.18) 

(iii) The equality Ve = Se holds on Ai. 

The proof of the last item of the previous theorem is achieved by characterizing 
the optimal path for the variational problem (|2.10p defining the quasi-potential. 
For E < Eq it is obtained by the following algorithm. Given p G Al let $(/9) € J-e 
be the solution to (|2l7l) and define G = e*(''V[l + e*^''^]. Let F = Ft{u) be 
the solution to the viscous Burgers equation (j2.5p with initial condition G and set 
■0 = log[i^/(l-F)], note that ■00 = ^{p) and 0* ^ as t ^ oo. Let p*^ = <^~'^{tpt), 
i.e. pI is given by the l.h.s. of (|2.17p with Lp replaced by ipt- Observe that p'g = p and 
Pj — > yo^ as t — > oo. The optimal path for (|2.10[1 is then ttj* = p*_f., the fact that it is 
defined on the time interval (— oo, 0] instead of [0, oo) makes no real difference. As 
discussed in [31 [5], this description of the optimal path tt* is related to the possibility 
of expressing the hydrodynamic limit for the process on E^r whose generator is the 
adjoint of ijv hi L2{dp^) in terms of (j2.5p via the nonlocal map $. 

The asymmetric limit. Consider the boundary driven asymmetric exclusion pro- 
cess, that is the process on S^r with generator given by (j2.ip where the external 
field E is replaced by Na and the generator is speeded up by N instead of N'^. 
We consider only the case a < 0. According to the previous notation, denote by 
Pno ^^'^ unique invariant measure of the boundary driven asymmetric exclusion 
process with external field aN. In the hydrodynamic scaling limit, it is proved 
in [Tj that the empirical density converges to the unique entropy solution to the 
inviscid Burgers equation with BLN boundary conditions, namely (|2.5p with E/2 
replaced by sinh(a/2) and no viscosity. 

Let p^ S {p-, p+, 1/2} be such that maxrg[p_ p^j x(^) = xiPa)- If is not difficult 
to check that the stationary profile p^ converges, as E — oo, to the constant 
density profile equal to p^, which is the unique stationary solution to the inviscid 
Burgers equation with the prescribed boundary conditions. 

By using combinatorial techniques, it is shown in [T3] that the sequence of prob- 
ability measures {p^^ o (tt^)"^} on A4 satisfies a large deviation principle with 
speed N and rate function S'a defined as follows. Let 

:= e Ci(hl, 1]) : (^(±1) = O} . (2.19) 

Note that J^e C J- a. 

Given p G M and ip G set 

Gaip.ip):^ y ^d7.{plogp+(l-p)log(l-p) + (l-p)<^-log(l + e^)-Aa} (2.20) 

in which the constant Aa is 

Aa := max logx(r) = logx(Pa) • (2-21) 

Let 

Saip) := sup Gaip,^) ■ (2.22) 
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t : gconv 



t : cgconv 



The functional Sa. is written in a somewhat different form in |14j . The above 
expression is however simply obtained by replacing the trial function F in [14j by 
e'^/(l + e"^). The advantage of the above formulation is that for each p G Ai the 
functional Ga{p,-) is concave on !Fa- By choosing (p = log[p^/{l — p^)] as trial 
function in (|2.22[) we get a lower bound analogous to (|2.15p : 

/I , 
du{plog^ + il-p)log-—^}. 

Note finally that Sa does not depend on a < 0. 

We prove in Section O that the functional Se converges, as i? | — oo, to 5'a. 
As discussed in [6l Lemma 4.3], the appropriate notion of variational convergence 
for rate functionals is the so-called F-convergence. Referring e.g. to [1^ for more 
details, we just recall its definition. Let X be a metric space. A sequence of 
functionals F„ : X [0, +oo] is said to T-converge to a functional F : X [0, +oo] 
if the following two conditions hold for each x ^ X. There exists a sequence Xn — > x 
such that lim„ Fn[xn) < (T-limsup inequality) and for any sequence Xn x 

we have lim„ i^„(x„) > F{x) {T-liminf inequality) . 

Theorem 2.4. Let Se ■ M ^ [0, +oo] be as defined in p.l4p . As E [ -oo, the 
sequence of functionals {Se} T-converges in M to S^ defined in (|2.22p . 

While the above result deals only with the variational convergence of the quasi- 
potential, it is reasonable to expect also the convergence of the dynamical rate 
functional. More precisely, the dynamic rate functional (j2.7p of the weakly asym- 
metric exclusion process should converge, in the appropriate scaling, to the one 
for the asymmetric exclusion process. We refer to [8] for a discussion of this topic 
and we mention that the above result has been proven in for general scalar 
conservation laws on the real line. 

F-convcrgcncc implies an upper bound for the infimum over open sets and a 
lower bound for the infimum over compacts sets: For each compact set IC d A4 and 
each open set C 



lim 



E- 



inf Se{p) 



> 



inf Sa{p) 



lim inf Se(p) 

E-,-oo peO 



< 



inf 5a (p) 

peo 



The proof of this statement is straightforward and can be found in 
Since A4 is compact, the previous fact and Theorems 12.21 l2T3] (iii). 
following asymptotics for the invariant measure p^. 



[H Prop. 1.18]. 
12.41 provide the 



Corollary 2.5. For each closed set C <Z Ai and each open set O d A4 



lim 

E^-oo 



1 



N IN 



EC) 



< 



- inf 5'a (p) 

pec 



lim 



lim T7logA'B 



N (N 



O) 



> 



- inf 5a (p) 
peo 



The last topic we discuss is the asymptotic behavior as, E —oo, of the solution 
to the Euler-Lagrange equation (|2.17p . More precisely, we show that it converges 
to the unique maximizer for (|2.22p . 

Consider the set J-a equipped with the topology inherited from the weak con- 
vergence of measures on [—1, 1): ^ in Ja if and only if J^^dif^ G ^ J^-^^dipG 



10 



L. BERTINI, D. GABRIELLI, AND C. LANDIM 



for any function G in Co([— 1, 1)), the set of continuous functions G : [—1, 1) ^ M 
such that hm„|i G{u) = 0. The closure of J-a, denoted by JFa, consists of all 
nondecreasing, cadlag functions if : [—1,1) — > such that 1) = </3-, 

lim„|i ip{u) < tp^. By Helly theorem is a compact Polish space. Moreover, if 
ip" ip in J-s, then ip"(u) -— > ^(u) Lebesgue a.e. 

t:cof Theorem 2.6. Fix p ^ Ad. There exists a unique (f> G .Fa such that >S'a(p) = 
maxipgjr^ C?a(p, = Qa{p,4>)- Let 4>e '■= '^{p) G &e the optimal profile for 
(j2.14p . j4s ^ — oo t/ie sequence {</)£;} converges to 4> in T^- 



s:5 



3. The nonequilibrium free energy 

In this section we analyze the variational problem (j2.14p and prove items (i) and 
(ii) in Theorem l2.3l We start by proving an existence and uniqueness result for the 
Euler-Lagrange equation (|2.17|) together with a dependence of the solution with 
respect to p. We consider the space C^([— 1,1]) endowed with the norm ||/||ci ■= 
ll/lloo + ll/'lloo where ||.g||oo := sup„g[_i \g{u)\. For each E < Eq the set 
defined in ()2.1ip is a convex subset of C^([— 1,1]); we denote by Te = {f E 
CHhl, 1]) : <^(±1) = <^± , > V £;} its closure in C^{[-1, 1]). 

t:deq Theorem 3.1. Let E < Eq. For each p Cz Ai there exists in Te o. unique solution 
to ()2.17p . denoted by $(jO). Furthermore, 

(i) //pgC([-1,1];[0,1]) then^p)eC^{[-l,l]). 

(ii) Let {p"} C A4 be a sequence converging to p in A4. Then {$(p")} C J^e 
converges to ^{p) in C^([— 1, 1]). 

Proof. The proof is divided in several steps. 

Existence of solutions. For E < 0, resp. E e {0,Eo), we formulate (|2.17p as an 
integral-differential equation informally obtained multiplying ()2.17p by ip' — E, resp. 
by ip' , and integrating the resulting equation. Existence of solutions will be deduced 
from Schauder fixed point theorem. 

Given E < Eq, p G M, and ip e J^e, let 

^'''(P,^;") [p-^-^][^'{u)~E], 

n^^\p,^-u) [p--_1-^](^'(k) . 

For a fixed p G and i 1,2 we define the integral-differential operators Kp : 
TE^C\[~l,l\)hy 



dvexpi^j dwTiP^\p,ip;w) 
J dv exp ^^J dwTZ'^^\p,ip]w) 



K.f\ip){u) := ^^+E{u + l) 



+ W+ 2E) 



J dv exp < 


j dwn^^\p,^;w)^ 




J dwTZ^^\p,^;w)^ 


J dv exp < 
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For E < Q, resp. E G (0, Eq), wc formulate the boundary problem (|2.17p as a fixed 
point on for the operator lC^p\ resp. /Cp^"*. 

Consider first the case E < corresponding to i = 1. Simple computations show 
that for each p G M the map ]C^p^ is a continuous on Te and IC'jP(^J-E) C J-e- It 
is also straightforward to check that there exists a constant Ci — Ci{(p-, (p+,E) e 
(0, oo) such that for any p G A4, ip G J-e, and u, v E [—1, 1], 



1 

Ci 



< 



d 



du 



dv 



du 



< Ci\u~v\. 

(3.1) 

In particular IC^^^^^J-e) C Te- Notice that Te is a closed convex subset of 
C^([— 1, 1]) and, by the previous bounds and Ascoli-Arzela theorem, /Cp^^ (-^-e) has 
compact closure in C^([— 1,1]). By Schauder fixed point theorem wc get that for 
each p E M. there exists ip* G Te such that K^p^(ip*) = p* . From (|3.ip it follows 
that p* G Te and standard manipulations show that p* satisfies (|2.17|) Lebcsgue 
a.e. 

The case E G (0, i^o), corresponding to a fixed point for /Cp , is analyzed in 
the same way. In this case, it is indeed straightforward to check that there exists 
a constant Ci = C2 (iy9_ ,(/?+, i?) G (0, 00) such that for any p G M., pi G Te-, and 
^G[-l,l], 



du 



-E<C2, ^/Ci^H(^)(w)-^/Ci"^M(ii) <C2\u-v\. 

(3.2) 

Uniqueness of solutions. Let (f> G Te, ^ 0, be a solution to (|2.17|) : by chain rule 
the equation 

1 , 6' -El' 6" 1 



E 



log- 



b' -E) 



= P 



holds Lebesgue a.e. Hence, for each u E [—1,1], 



E 



log 



(f)'{u) - E 
0'(w) 



E 



log 



4''i-l)-E 
0'(-l) 



dv 



p{v) 



1 



(3.3) 



If 4)'i{ — l) — '/'2(^l) an application 
We next assume 0i(— 1) < 02 (—1) 



Let (^^i , 02 G -^B be two solutions to (|2.17p . 
of Gronwall inequality in (|3.3p yields 0i ~ 02. 
and deduce a contradiction. Recall that 0^ > V and let u := inf{w G ( — 1, 1] : 
(j)i{v) — 02(i^)} which belongs to (—1,1] because 0i(±l) = 02(±1) and 0'j(— 1) < 
02 (—1). By definition of u, 0i(u) < 02 (m) for any u G ( — l,u), 0i('ii) = 02 (u) and 
0i(M) > 02(")- Note that the real function (0 V 00) 9 z ^ E^^ ^og[{z - E)/z\ 
is strictly increasing. Therefore from (|3.3p we obtain (p'liu) < 02 (m), which is a 
contradiction and concludes the proof of the uniqueness. 

The case E = can be treated similarly, with —(1/0')' in place of {{1/E) log[(0' — 
E) /((>']}', and was examined in [J. 

Claims (i) and (ii). Claim (i) follows straightforwardly from the previous anal- 
ysis. To prove (ii), let 0" $(p") G Te- By and Ascoli-Arzela 
theorem, the sequence {0"} C Te is precompact in C^([— 1,1]). It remains to 
show uniqueness of its limit points. Consider a subsequence rij and assume that 
{(/)"'} converges to ijj in C^([— 1,1]). Since {p"'} converges to p in and {0"^} 
converges to tp in C^([— 1,1]), for E < 0, resp. for E G (0, i?o), we have that 
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/C^ij (0"^ ) converges to ICp^\ip), resp. IC^\}j {(t>"'^) converges to icf\'ip). In particu- 
lar, ip = limj = limj /c|,nj (0"^ ) = ^p'lV') for i = 1, 2. By the uniqueness result, 
ip ~ ^{p). This shows that $(jo) is the unique possible limit point of the sequence 
{((>"}, and concludes the proof of Claim (ii). □ 

Fix a path p EE pt{u) G C^'"{[0,T] x [-1, 1]; [0, 1]) and let <j) ee ^{pt){u) be the 
solution to dSUl). By TheoremO belongs to C^-'^{[0,T] x [-1,1]). Note also 
that, by (|3.1[) and (|3.2p . for each E < Eq there exists a constant C G (0,c») such 
that for any {t,u) G [0,T] x [-1, 1] 

C^' < VMu) < C ifE<0, , 

< V(j)t{u)-E < C ifO <E<Eo. ^ ^ 



t:ddeq| Lemma 3.2. Let E < Eq, T > 0, p £ C^^°{[0,T] x [~l,l];[0,l]), and (j) := ^{pt) 
be the solution to (|2.17p . Then (p G C^'^([0,T] x [—1, 1]) and ip := dtcj) is the unique 
classical solution to the linear boundary value problem 



^Pt=dtpt G [0,r] x (-1,1) 



(l + e-^*)^ ' ' (3.5) I dPeqf 



-V0t(V0t -S) 
^t(±l) = te[0,T] 



Proof. Fix t G [0, T]. For h ^ such that t + h £ [0, T] define Vt''(-) by ip^iu) := 
[4>t+hiu) - Mu)]/h. By Theorem 0(1), i'ti-) belongs to C2([-l,l]). Set R'l := 
[pt+h — Pt\/h; from (|2.17p it follows that tp'^ solves 

Ay.j' A(j)t+h{y(pt + ^(pt+h-E) 

- E) V0,(V0, - i?)V0,+,(V0,+, - E) ^^-^^ 



p0t p'l'At _ 1 



(l + e'^*)(l + e'^*+'') h 

for (t,M) G [0,r] X (-1, 1) with the boundary conditions Vt''(±l) = 0, t G [0,T]. 

Multiplying the above equation by ip^ and integrating in du, using the inequality 
x(e^ — 1) > and an integration by parts we get that 

VcPt{V^t-E), 

where 



F{(ptAt+h) 



pel 



(V0O'(V</)t - EfWcPt+hi'^cPt+h - E) 
X {A(/)tV0t+,,(V0t+ft - ^)(2V0t - E) 

- A(Pt+hyM^<Pt - S)(V(/)t+,, + V(/)t - S)} . 

For each t G [0,T], 

lim ||^^(0t,0t+ft)|loo = . (3.8) [f^ 

/i^O 

Indeed, since p G Ci'O([0,T] x [-1,1]), as /i ^ 0, pt+/i(-) ^ /9t(-) in C([-l,l]). 
By Theorem 13.11 fii). (pt+h{-) — * 0t(') in C'"'^([~lj !])• By the differential equation 
(fSTTT]) . 0t+/i(-) ^ 0t(-) in C2([-l, 1]). Together with ((3^ this concludes the proof 

of dsT 
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By p.4p . Cauchy-Schwarz, and Poincare inequality for the Dirichlet Laplacian 
in [—1,1], we obtain from (j3.7|) that 

< c'(vv'^v^f)l/2 [{RtR'ty/' + \\F{ci,t,ci,t+,,)\\^{vit^^t)'^' 

for some constant C > 0. 

From (13.91) and (13.81) it follows that there exists a constant C" > such that 



csp 



lim {V^j^, Vij^) < C" {dtpt,dtpt), te[0,T]. (3.10) ubhl 

h^O 

Therefore for each t G [0,T] the sequence {ipti')} precompact in C([— 1, 1]). By 
taking the limit h in (j3.6p and using (|3.8p . it is now easy to show that any 
limit point of {i/Jti-)} is a weak solution to p.Sp . By the classical theory on one- 
dimensional elliptic problems, see e.g. [211 IV, §2.1], there exists a unique weak 
solution to (|3.5p which is in fact the classical solution because dtpt{-) belongs to 
C([— 1, 1]). This implies that there exists a unique limit point '0t(-) G C^([— 1, 1]). 
Finally £ C*''^([0, T] x [-1,1]) by the continuous dependence in the C^([-l,l]) 
topology of the solution to (|3.5p w.r.t. dtpt{-) in the C([— 1, 1]) topology. □ 

We are now in a position to prove two statements of the first main result of this 
article. 

Proof of Theorem \2.3\ (i) and (ii). We start with Claim (i). The case E — Eq 
follows from the definition (|2.16p of the functional Sec- Assume E < Eq. By the 
convexity of the map xip) ~ plogp + (1 — p) log(l — p), for each ip e J-e the 
functional Ge{-,^) is convex and lower semicontinuous on A4. Hence, by (j2.14p . 
the functional Se, being the supremum of convex lower semicontinuous functionals, 
is a convex lower semicontinuous functional on A4. On the other hand, since the 
real function (0V£^,oo) 3 x \_x log x — {x — E) log(.T ~ E)\/ E is strictly concave, 
Jensen inequality and <p(±l) = ^p± imply that Ge{p, is bounded by some constant 
depending only on ip± and E. This proves (i). 

Fix p E M. The strict concavity mentioned above and the strict concavity of 
the real function M 9 x t-^ — log (l -|-e^) yield that the functional Ge{p, •) is strictly 
concave on J-e- Thanks to Theorem 13.11 it easily follows that the supremum on 
the r.h.s. of (|2.14p is uniquely attained when ip = $(jO). □ 

In the proof of the equality between the quasi-potential Ve and the functional 
Se, we shall need the following simple observation. 

t:Sden Lemma 3.3. For each p Cz M there exists a sequence {p"} C A4 converging to p 
inM and such that: p" £ C2([-l,l]), p"(±l) = p±, < p" < 1, S'£;(p") ^ Se{p). 

Proof. For E ~ Eq, this is obvious from the definition of the functional Seo- For 
E < Eo, given p e M, it is enough to consider a sequence {p"} C C^([— 1,1]) 
with p"(±l) = p± and < p" < 1, which converges to p du a.e. By Theorem 12.31 
(ii), Theorem 13.11 fii). and dominated convergence, Se{p'^) = ^_e(p", $(p")) — > 

gE{p,Hp)) ^Se{p). □ 
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4. The quasi-potential 



In this section we characterize the optimal path for the variational problem (|2.10p 
defining the quasi-potential Ve and conclude the proof of Theorem 12 . 3 1 by showing 
the equality Ve ~ Se- The heuristic argument is quite simple. To the variational 
problem ()2.10|) is associated the following Hamilton- Jacobi equation [3l [5]. The 
quasi-potential Ve is the maximal solution to 

'SVe 1. E_ , „ ^^^^^ 



1 / SVe ^yy^^EX ^ /OVe i . 







with the boundary condition that 5Ve/5p vanishes at the endpoints of [—1, 1]. Few 
formal computations show that Se solves (j4.ip . To check that Se is the maximal 
solution one constructs a suitable path for the variational problem (|2.10p . [H [5] . 
Since it is not clear how to analyze (|4.ip directly, we first approximate, as in [1], 
paths TT G D([0,T]; A^) with It{'t^\Pe) < oo by smooth paths bounded away from 
and 1 which satisfy the boundary conditions p± at the endpoints of [—1,1]. For 
such smooth paths we can make sense of ()4.ip and complete the proof. 

In the case E = Eg, the process is reversible and the picture is well known. 
The path which minimizes the variational formula defining the quasi-potential is 
the solution of the hydrodynamic equation reversed in time. The identity between 
Seo Emd Veq follows easily from this principle. The proof presented below for 
E < Eq can be adapted with several simplifications. It is enough to set $(p) = 

loglPiJo/l - Peo} everywhere. 

Assume from now on that E < Eq. We first need to recall some notation 
introduced in [7]. Fix a density profile 7 : [—1, 1] — > [0, 1] and a time T > 0. Denote 
by J^2 = ^2{T, 7, p±) the set of trajectories n in C([0, T],M) bounded away from 
and 1 in the sense that for each t > 0, there exists e > such that £<7r<l— eon 
[t,T]; which satisfy the boundary conditions, ttq = 7, 7rt(±l) — /0±, < t < T; and 
for which there exists 61, S2 > such that ttj follows the hydrodynamic equation 
(j2.5p in the time interval [0, Si], TTt is constant in the time interval [61, 61 + 62] and 
TTt is smooth in time in the time interval ((5i, T]. 

If the density profile 7 is the stationary profile p^, the trajectories tt in J- 2 are in 
fact constant in the time interval [0, Si + ^2]. Since they are also smooth in time in 
{5i, T], we deduce that they are smooth in time in the all interval [0, T]. Moreover, 
since Pe is bounded away from and 1 , there exists e > such that e < tt < 1 — e 

on [o,r]. 

Assume that J = Pe and recall from the proof of [71 Theorem 4.6] the definition 
of the sequence of trajectories {tt^ : e > 0}. Since a path tt in .F2 is in fact constant 
in the time interval [0, 6], each tt^ is smooth in space and time. In particular, let 



._ (700,00 ^jQ^y] X [-1,1]) r\T2 ■ 

Theorem 4.6 in [7j can be rephrased in the present context as 



(4.2) 



Theorem 4.1. For each tt in D{[0,T], M) such that It{tt\pe) < oO; there exists a 
sequence {tt"} C T>o converging to tt in I?([0, T]; yVf) such that ^T(7r"|p^) converges 
to It{tt\pe)- 

The first two lemmata of this section state that, for smooth paths, the functional 
Se satisfies (|4.ip . Recall that for p £ Al we denote by ^(p) e !Fe the unique 
solution to (|TT7)1 . 



ham-jacO 



dDd 
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dots I Lemma 4.2. Let E < Eq, T > 0, n e Vq, and T : [0,T] x [-1, 1] ^ R 6e defined 
by 

Ft := \og—^ $(7rf) . (4.3) I dGt 

1 - TTt 

Then 



Se{ttt) - Se{tto) = I dt{Tt,dtTTt). (4.4) | nome2 

JQ 



Proof. Let (j) EE (j)tiu) $(7rt) (w), {t,u) G [0,T] x [-1,1]. By Lemma | 
belongs to Ci'2([0,r] x [-1,1]). Since 0t(±l) = (p±, then 9t0t(±l) ^ 0, t e 
[0,T]. By Theorem 12.31 (ii). dominated convergence, an exphcit computation, and 
an integration by parts, 

P.in.) = I ^.(.., $(..)) 

The lemma follows noticing that the last term vanishes by (|2.17p . □ 
Let 

TUo := {/5GC2([-1,1]) : p(±l)=p±, 0<p<l} . (4.5) [dM 

ham-jac Lemma 4.3. Let ii^ < i?o, p G A^O; a?^'^ T : [—1, 1] ^ M. be defined by 

r := log— $(/?) . (4.6) [hTo" 

1 - p 

Then, 

(vr , x(p) vr) - (Vp - i?x(p) , vr) = o . (4.7) [IT 

Proo/. As before we let (/> = 0(u) := <I>(p) (u), u G [-1, 1]. By Theorem [3T] (i) . (j) 
belongs to C^([— 1, 1]). By the definition of F in (|4.6|) . statement ()4.7|) is equivalent 
to 

(Vp, -V0 + i?) + (-V(/., x(p)(-V0 + i?)) = 0. 
The above equation holds if and only if 

^(P-T^) ' ^^-e) + (^(t^) ' ^^-e) - (^'^' x(p)(V0-i?)) = . 



l + e"^/' / \ Vl + e"^. 

Since e-^f^iVfl + e"^*^^^] = e'^*/[l + e'^±] = p± = /o(±l), integrating by parts the 
previous equation, it becomes 

p-j^,A^)- {{j^ - x(p))v0,V<^-i.) ^ 0. 

At this point the explicit expression for x given by xip) = /o(l ^ p) pla-ys a crucial 
role. Indeed, for such x, 

so that (14. 7p is equivalent to 

(p-^,A0 + V0(V0-i.)(l-p-^)) = 0, 
which holds true because = $(p) solves p.l7p . □ 
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We next prove the first half of the equahty Ve = Se- In fact the argument 
basicaUy shows that any sohition to the Hamihon-Jacobi equation ()4.1|) gives a 
lower bound on the quasi-potential. 



Proof of Theorem \2.S[ the inequality Ve > Se- In view of the variational definition 
of Ve in (|2.10[) , to prove the lemma we need to show that for each p G Ai wc have 
Se{p) < It{t^\Pe) foi' T > Q and any path tt e -D([0, T]]M) such that ttt = p. 

Assume firstly that p G A^o and consider only paths tt £ Dq- Of course the 
energy Q(7r) of such a path tt is finite. In view of the variational definition of 
It{t^\Pe) given in (|2.6p . (|2.7p . to prove that Se{p) < It{t^\Pe) is enough to exhibit 
some function H e Cl'^{[0,T] x [-1, 1]) for which Se{p) < JT,H,p^{Tr). We claim 
that r given in (|4.3p fulfills these conditions. Let 4> = (f>t{u) ■= ^{T^t) (u). Since tt S 
Vo, by LemmalOr € C^-^{[0,T] x [-1, 1]). On the other hand, since 7rt(±l) = p± 
and (/)t(±l) (^±, rt(±l) = 0, i e [0,T]; whence T S Co'^([0,T] x [-1,1]). 
Recalling the definition of the functional Jt.t.p^, after an integration by parts, we 
obtain that 

Jtt.p-Jtt) = ^ dt[{Tt,dt7rt) + ^{VTt,S77rt~Exi7rt))~^{xi7Tt),i^'rty) 

From Lemmata 14 . 2 [ and since ttq = p^, Se{Pe) = 0) it follows that JT,r,p£.(7r) ~ 
Se{p), which proves the statement for p G Mq and paths tt £ Pq. 

Let now p G A4 and consider an arbitrary path tt G D(\0,T]; Ai^ such that 
ttt = p- With no loss of generality we can assume JT(7r|/0^) < oo. Let {tt"} C Vq 
be the sequence given by Theorem 14.11 The result for p S A^o and paths in Vq, 
together with the lower semicontinuity of Se, yield 

lT{Tr\pE) = lim It{t^'^\Pe) > liffi Se{tit) > Se{t^t) = Se{p) , 
which concludes the proof. □ 

To prove the converse inequality Ve < Se on A^, we need to characterize the 
optimal path for the variational problem (j2.10p . The following lemma explains 
which is the right candidate. 

Denote by C'^{Q.t) the smooth functions H : Q.t R with compact support. 
For a trajectory tt in D{[0,T],M), let TiKxi'^)) be the Hilbert space induced by 
C^(J7r) endowed with the scalar product defined by 

((G,i/))i,^(,) ^ dt J ^du{\/G){t,u){\/H){t,u)x{At,u)) . 

Induced means that we first declare two functions F, G in C^(57t) to be equivalent 
if {{F — G,F — G)}! ,^(7r) = and then we complete the quotient space with respect 
to scalar product. Denote by || • ||i.T(.(7r) the norm associated to the scalar product 

Repeating the arguments of the proof of Lemma 4.7 in [7], we obtain an explicit 
expression of the rate function /T(vr|7) in terms of a solution to an elliptic equation. 

Lemma 4.4. Fix a trajectory tt in Vq- For each < t < T , let Ht be the unique 
solution to the elliptic equation 

r dtTTt = (1/2) A^t ~~ V{x{7Tt)[{E/2) + \/Ht]} , , 
\ Ht{±l) = . ^^-^^ 
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Then, H is smooth on [0,T] x [—1,1] and 

lT{7T\no) = ImiUi.) • (4.9) ^ 

We could have used next lemma to prove the inequality Ve > Se', we presented 
the separate argument before for its simplicity. On the other hand, (|4.11[) clearly 
suggests that the optimal path for the variational problem (j2.10|) is obtained by 
taking a path which satisfies ()4.10p with K = 0. Recall that <i>(p) denotes the 
solution to ((TT7)l . 

I=I*+S Lemma 4.5. Let E < Eq, T > 0, 7 G Mq, and t: e Vq such that /T(7r|7) < 00. 

Then, there exists K in Cq'^([0, T] x [—1,1]) such that n is a classical solution to 

'dtTTt + I Vx(7rt) = -i ATTt + V[x(^t)V($(^t) + Kt)] 

7rt(±l) = p± (4-10) I f 07bis 



.TTO = 7 

Furthermore, 



lTi7rh) = SEi7:T)~SEij) + l\\K\\l ■ (4.11) ^ 



2 II ''i,xM 

Proof. Note that 7 = ttq because we assume the rate function to be finite. Denote 
by H the smooth function introduced in Lemma 14.41 and let F be as defined in 
()4.3p . We claim that K := T — H meets the requirements in the lemma. As before 
we have that F belongs to Co'^([0,r] x [—1,1]). Hence, K also belongs to this 
space because H is smooth and vanishes at the boundary of [—1,1]. The equation 
(|4.10p follows easily from replacing H hy T - K. To prove identity (|4.1ip . 

consider (14. 4p and express dtiTt in terms of the differential equation in (|4.8p . Since 
H = r — K , after an integration by parts we get that Se^t^t) — SEi'j) is equal to 

-i / dt{\/rt,\/nt^Ex{7rt)) + I dt{WTt,x{T^t)y{Tt-Kt)) . 
^ Jq Jo 

By Lemma 14.31 the previous expression is equal to 

di(VFt, x(7rt)VFt) - / dt{VTt,x{^t)S7Kt) . 
^ Jo Jo 

Since K = T — H, we finally get that 

1 2 1 ||2 

Se{7tt) - SEh) + 2ll^lli,xW = 2"^"i.xW 



which, in view of (|4.9p . concludes the proof. □ 

We next show how a solution to the (nonlocal) equation (|4.10p with A' = 
can be obtained by the algorithm presented below the statement of Theorem 12.31 
Recall that such algorithm requires to solve (|2.17p only for the initial datum and 
then to solve the (local) hydrodynamic equation (|2.5p . Note indeed that by setting 
TTj := p!L(, where p* is defined in the next lemma, then tt* solves the differential 
equation in (|4.10|1 with K ^0. 

Fix E < Eo, e Mo and set G e*(T)/[l + e*(T)]. By Theorem O the profile 
G belongs to C^([— 1, 1]), it is strictly increasing and satisfies G(±l) = p±. Denote 
by F = Ft{u) G C^'^([0, 00) x [—1, 1]) the solution to the hydrodynamic equation 
(|2.5p with 7 replaced by G. By the maximum principle, p_ < F < 
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t : cadh I Lemma 4.6. Let -0 := log[F/(l - F)]. Then, belongs to C^'^{[0, oo) x [-1, 1]) 
and satisfies Vij} > V i?. Let p* = p1{u) be defined by 

Then, p* belongs to C^'^([0, oo) x [—1, 1]), satisfies (±1) = p±, < p* < 1, and 
solves 

'dtp; ~ f vx(p:) - 1 Ap; y[x{p:)vHp:)] 

p;{±l)=p± (4-13) Udhr* 

.Po = 7 ■ 

Proof. Let t/j : [0, oo) x [—1, 1] ^ M be given by ip = log{F/l — F} and set 

T := sup{t > : VV's(it) > V S for all {s,u) G [0,i] x [-1,1]} . 

Since VV'o = V<i>(7) > V t > by continuity. We show at the end of the proof 
that r = oo. 

A straightforward computation shows that ijj solves 

1 1 - 




VV'(V?/; - E) 

(4.14) 



Since V e C'^''^{[0,oo) x [-1,1]), definition KWi yields p* € C^-^{[0,t) x [-1,1]) 
and Pq ~ 7. On the other hand, from (|4.14p we deduce that for any t G [0, r) 

1 — e*^* 

A^t(±l) + ^__V7^t(±l) [VV;t(±l) - £;] = . 
Whence, again by (|4.12p . 

1 1 - e"^* 

1 + e'^i 1 + e'^i 

By using (|4.14p . a long and tedious computation that we omit shows that p^, 
t G [0, t), solves the differential equation in (|4.13p . 

We next show that < p* < 1. Since 7 G A^o, there exists S G (0, 1) such that 
5 < ^ < 1 — 6. We claim that mm{p^, 1 — p+,S} < p* < ma.x{p^, 1 — p-,1 — 6}. 
Fix t G (0, r) and assume that (•) has a local maximum at uq G (—1, 1)- Since p* 
solves (jLTSl) . since $(/3*) solves ()2.17p and since Vp^ (mq) = 0, Ap^^uo) < 0, 

dtpUuo) = ^Ap;iu„) - xip*tM)A'l>ip;){uo) 

< -x{p*t{uo)) V$(p:)(uo) [V$(p:)(^.o) - E] [ptiuo) - ^ ^ J^,,)^^„y 

Assume now that Pt (mo) > 1 — P-- Since ^{pt ) > (f- we deduce (wo) — [1 + 
g<E.(p:)(uo)]-i > l-p_-[l + e'^-]-i = 0. As V$(pt*)(uo) > OV^;, we get atPt*(uo) < 
0. In particular, by a standard argument, p* < max{p_|-, 1 — p_, 1 — 6}. The proof 
of the lower bound is analogous. 

We conclude the proof showing that r = 00. Assume that r < 00. Since for each 
t G [0,r), p* belongs to Mo, it follows from that $(p*) = Vt, t G [0,t). By 

TheoremlOfii). $(p;) = i/v so that Vi/)^ = V$(p;) > i^VO because $(p*) belongs 
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to Te- By continuity, there exists (5 > such that Vi/'t > i? V for r < t < r + (5. 
This contradicts the definition of r. □ 

Fix a density profile 7 : [—1, 1] [0, 1], a time T > and consider the solution 
p* to (|4.13p . Let A((-) = PT-t(')- Clearly, A is the solution to (|4.10|) in the time 
interval [0,T] with if = and initial condition Ag = Px- In particular, by (|4.1ip . 

It{Mpt) = Se{i) - Se{pt) ■ 
In the next lemma wc prove that converges to as T ^ cxd. Letting T | 00 
in the previous formula, we see that the time reversed trajectory of (|4.13p is the 
natural candidate to solve the variational formula defining the quasi-potential. This 
argument is made rigorous in the next paragraphs. 

By standard properties of parabolic equations on a bounded interval, see e.g. 
[TT] . as t ^ 00, the solution to (|2.5p converges, in a strong topology, to the unique 
stationary solution p^. Such convergence implies that the path p* , as defined in 
Lemma l4!6l also converges to p^; as t — > 00. This is the content of the next lemma. 
This result will permit to use the time reversal of p* as a trial path in the variational 
problem ((2TT0)) . 

t:conv Lemma 4.7. Let E < Eq, 7 e -Ado, o,i^d p* be defined as in Lemma \4-(>\ A.s 
t — > cxo, the profile pi G A^o converges to p^ in the C^([— 1, 1]) topology, uniformly 
for-f e Mq. 



Proof. Recall the notation introduced just before Theorem l4.6l Let p be the solution 
to (|2.5p . In [Til Theorem 4.9] it is shown that, as t ~* 00, the profile pt converges 
to in the C"'^([— 1, 1]) topology, uniformly for 7 G Aio- By the methods there 
developed, it is however straightforward to prove this statement in the C'^([— 1, 1]) 
topology. In particular. Ft converges to p^ in the C'^([— 1,1]) topology to p^ so 
that V't converges to log[p£/(l + p^)] = (p^ in the C'^([— 1, 1]) topology uniformly 
in 7 G A^o- Since ^{pe) = 'pE^ the statement now follows from (|4.12p . □ 

We next show that profiles close to p^ in a strong topology can be reached with 
a small cost. 

t:join Lemma 4.8. Let E < Eq andS G (0, 1). Then, there existT > and constantC = 
C{E, p±, 6) > such that the following hold. For each p G C^([— 1,1]) satisfying 
p(±l) = p± and S < p < 1 — S, there exists a path n G D(\0,T]; A4) such that 
ttt = p and 

II 1 1 2 

/t(7t|pb) < C\\p-Pe\\^, . 

Proof. Simple computations show that T = 1 and the straight path ttj = p^ + t{p — 
Pg) meet the requirements. For i? = 0, in [U Lemma 5.7] a more clever path is 
chosen which yields a bound in terms of the L2 norm of p — p^. □ 

We can now conclude the proof of Theorem [ 



Proof of Theorem \2.S[ the inequality Ve < Se- Given p G and (5 > we need 
to find T > and a path n* G D{[0,T];M) such that tt^ = p and It{t^*\pe) < 
Se{p) + S. By Lemma 1531 there exists a sequence {p"} C A^o converging to p in 
M and such that Se{p"') — > Se{p)- Let p*'" be the path constructed in Lemma l46l 
with 7 replaced by p" and pick e > to be chosen later. By Lemma 14771 there exists 
a time Ti = Ti(e) > independent of n such that ||p^" — PeWc^ < £■ Whence, 
by Lemma 14751 there exists a time T2 > 0, still independent of n, and a path tt", 
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t G [0,T2] such that ttq = Pe- ~ PtI^ ^^"^ ^T2('^"Ip_e) ^ Pe, where (3^ vanishes 
as £ — > and is independent of n. We now set T := Ti + T2 and let tt^ '", t G [0, T] 
be the path defined by 

*,n ^ /^rr i G [0,T2] 

which satisfies tTq'" = and tt^'" = p". The covariance of / w.r.t. time shifts, 
Lemmata 14.51 and 14.61 yield 

It[^*-'^\Pe) = ItA*^\Pe)+ItAPt:-\Pt:) 

< Pe + Ssipn - Se{Pt;') < f3e + Se{p") ■ (4-15) ^ 

Since SeIp"^) Se{p) < oo and It{'\Pe) compact level sets, see Theorem 12. 1[ 
the bound (|4.15p imphes precompactness of the sequence {tt*'"} C D(J0,T];M). 
It therefore exists a path tt* and a subsequence nj such that tt*'"^ — > tt* in 
D([0,T];A^). In particular tt^ = linij tt^'"^ = limj p"^ = p. The lower semi- 
continuity of It{'\Pe) ^i^d (|4.15p now yield 

It{7t*\pe) < lim It{7t*'''^\pe) < Pe + Hm 5b(p"0 = /3s + ^^(p) 



which, by choosing e so that /?£ < (5, concludes the proof. □ 



s : asyl 



5. The asymmetric limit 

In this section we discuss the asymmetric limit E ^ ~oo and prove Theorems 
and 



Proof of Theorem \2.4\ T-liminf inequality. Fix p G M and a sequence {pe} C A4 
converging to p in as £' — cxo. We need to prove that lim ^ Se{pe) > Si^^p). 
Let Je be such that (|2.8[) holds; it is straightforward to check that 

lim — ^ = max x{^) j 

S^-oo re[p_,p+] 

whence, recalling that has been defined in (|2.13p and A^, in (|2.2ip . 

lim [Ae -log{-E)] = max logxW = . (5-1) 

-E^-oo re[p_,p+] 

Fix (p e C^+^([-l, 1]) such that (p(±l) = (p± and ip' > 0. From ((57T|) it easily 
follows that 



lim 

E-t~o 



J du[^[^' log ^' -{^' ~E)log{^' -E)] -{Ae-A,)}=0. (5.2) 

Recalling (|2.14p . (|2.12p and (|2.20p . from the convexity of the real function F : 
[0, 1] ^ R, F(p) = plogp + (1 - p) log(l - p) and ([521) we get 

liin Se{pe) > lim Ge{pe,^) > Gi,{p,^) ■ 

The proof of the F-liminf inequality is now completed by optimizing on (p. Note 
indeed that the supremum in (j2.22p can be restricted to strictly increasing (p e 
Ci+i([-l, 1]) such that (p(±l) = (p±. □ 
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Proof of Theorem \2.4\ T-limsup inequality. Fix p G A^, we need to exhibit a se- 
quence {pe} C M converging to p va M as E ^ — oo such that \\mE Se{pe) < 
Sa{p). Wc claim that the constant sequence pE = p meets this condition. 

Recahing item (ii) in Theorem l2.3i let (j)E ■= ^{p) € ^e be the solution to (|2.17p 
in which we indicated explicitly its dependence on E. From the concavity of the 
real function F : [0, oo) R, F{x) = E''^ [xlogx - {x - E) \og{x - E)], E < 
Jensen inequality, and (|5.ip we deduce 

lim 

E^-o 



J " ^"^'^ " ^^S^*^'^ " - ("^E - A,)} < . (5.3) I pezzo2 



Since J-e C J-g, and is compact, the sequence {0b} is precompact in Let now 
0* G ^a be any limit point of {0b} and pick a subsequence E' —oo such that 
4>E' — * 0* in -^a- In particular 0b'(u) ^ 0*('u) Lebesgue a.e. Recalling Theorem l2.3l 
(ii), ((TT21) . (|2:20l) . and using jSS]) we get that 

lim Se'{p)^ lim Ge' {p, (I^e') < Ga{p, (/)*)< Si,{p) , 

E — » — oo — » — oo 

which concludes the proof. □ 

Proof of Theorem \2.6[ Existence of a maximizer for (|2.22[) follows from the com- 
pactness of and from the continuity of Ga{p, ■) for the topology of Ta. On 
the other hand, the strict concavity of the function F : R+, F{ifi) = 

— log(l -I- e'^), gives the uniqueness of the maximizer. 

The proof of the convergence of the maximizers follows a variational approach. 
Given p G M and < we define Ge{p, ■) ■ Ta ^ [— oo, +oo) by 



Ge{p-, <p) 



\Ge{p, ^) if G Te 
1 — cx) otherwise. 



By [ini Theorem 1.21], with all inequalities reversed since we focus on maximizers 
instead of minimizers, the convergence of the sequence {0b} to in Ta follows 
from the next three conditions. Fix p G Ai and Lp G Ta then: 

(i) for any sequence v^b ^ (ys in J?^a, liniB Ge{PtVe) < Gaip, ^); 

(ii) there exists a sequence ipE V' in -^a such that lim ^ Ge{Pt'Pe) > Ga.{p^ </?); 

(iii) is the unique maximizer for the functional Ga.{p, ■) on Ta- 

Proof of (i). We may assume that (/Sb G ^e] the proof of (i) is then achieved by 
noticing that (|5.3p holds also if 0b is replaced by v3b- 

Proof of (ii). Assume firstly that belongs to C^([— 1, 1]) and satisfies <i3(±l) = Lp±, 
kp' > 0. Since (|5.2p holds for such ip, it is enough to take the constant sequence ipE = 
(p. The proof of (ii) is completed by a density argument, see e.g. [TOl Rem. 1.29]. 
More precisely, it is enough to show that for each (p G Ta there exists a sequence 
V?" G Ci([-1, 1]) satisfying v?"(±l) = (^±, {ip")' > and such that ip"" ^ ip in Ta, 
GaiPi'P^) GaiPi'p)- This is implied by classical results on the approximation of 
BV functions by smooth ones. 

As we have already shown (iii), the proof is completed. □ 
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